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ABSTRACT 

We present an analysis of the Suzaku observations of the prototype wind- 
blown bubble NGC 6888 which is based both on use of standard spectral models 
and on a direct comparison of theoretical models with observations. The X-ray 
spectra of NGC 6888 are soft and most of the X-rays are in the (0.3 - 1.5 keV) 
energy range. But, hard X-rays (1.5 - 4.0 keV) are also detected (~ 10% of the 
observed flux). The corresponding spectral fits require a relatively cool plasma 
with kT < 0.5 keV but much hotter plasma with temperature kT > 2.0 keV 
is needed to match the observed hard X-ray emission. We find no appreciable 
temperature variations within the hot bubble in NGC 6888. The derived abun- 
dances (N, O, Ne) are consistent with those of the optical nebula. This indicates 
a common origin of the X-ray emitting gas and the outer cold shell: most of the 
X-ray plasma (having non-uniform spatial distribution: clumps) has flown into 
the hot bubble from the optical nebula. If the electron thermal conduction is ef- 
ficient, this can naturally explain the relatively low plasma temperature of most 
of the X-ray emitting plasma. Alternatively, the hot bubble in NGC 6888 will 
be adiabatic and the cold clumps are heated up to X-ray temperatures likely 
by energy exchange between the heavy particles (hot ions diffusing into the cold 
clumps). 
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Introduction 



Optical nebulosities trace large cavities around early-type stars (O, Of and Wolf-Rayet, 
WR) that are formed as the high-speed wind sweeps up ambient interstellar gas and com- 
presses it into a thin shell. The flow pattern, resulting from this interaction, was first 
recognized by Pikelner (1968) and consists of two regions of shocked gas: one corresponding 
to the shocked stellar wind and the other to shocked interstellar gas. The two regions are 
separated by a contact discontinuity and the hot interior gas can cool via thermal conduction 
across the interface. Because of its high temperature, the interior of the wind-blown bubbles 
(WBB) is expected to emit in X-rays. 

In 1970s and 1980s, analytical solutions for the WBB structure were derived that re- 
vealed details in the physics of these interesting objects. For a review of the analytical works 
and the related physical ideas see Dyson (1981) and McCray (1983). It should be noted that 
the most complete (semi-) analytical study of the WBB structure was presented in Weaver et 
al. (1977). With the increasing computing power, the numerical hydrodynamic simulations 
became a powerful tool for study i ng the physics of these objects. The first numerical mod- 



conductive WBB were studied numerically (e.g., 
1985a , b,c; Brighenti & D'Ercole, 1 995a,b, 1997; 



Rozvczka 


1985 


D'Erco 


c 


1992 



Rozvczka fc Tenorio-Tagle, 



Garcia- Segura fc Mac Low 



1996bl ; IZhekov fc Myasnikov 



1998 



19951 ; iGarcia-Segura et al.l Il996at iGarcia-Segura et al. 
2000). Thanks to numerical simulations, it was found that the WBB are subject to nu- 
merous dynamic instabilities. These simulations also allowed a much more complete and 
complex physical picture be explored, namely, by following various pha ses of the WBB evo- 



lution when the wind of the central stars varies with time (for details see lGarcia-Segura et al. 
1996at IGarcia-Segura et al.lll996bl ). 



Observational properties of the optical nebula in WBB were described in detail in the 
early works of Chu et al. (1983 and the references therein) and Lozinskaya (1982; see also 
the book by Lozinskaya 1993 for an observational review). But, it should be emphasized 
that the physics of the hot interior (hot bubble) of a wind-blown bubble is a cornerstone 
in the entire physical picture of these objects. This is so since the hot bubble is the region 
where the stellar wind energy is stored and subsequently used for driving the entire structure. 
From such a point of view, X-ray observations of WBB are very important because they can 
provide us with details about the physical conditions in the hot bubble. 

The first su ccessful X-ray detection of a WBB was that of NGC 6888 by Einstein 
( iBochkarevi Il988l ) . It was found that the X-ray emission from NGC 6888 is characterized 
by a plasma temperature kT = 0.28 — 0.8 keV (90% confidence interval), or T = 3.2 — 
9 MK. ROSAT observations were sensitive to cooler plasma and yielded a characteristic 
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temperature T pa 2 MK. They sh owed that the X-ray emission arises primarily in filament 
like structures ( IWrigge et al.lll994j ). ASCA observat ions suggested that even hotter plasma at 
T pa 8 MK might be present in the bubble interior (IWrigge et al.l 120051 ). Preliminary results 
from a recent Chandra observation, which covered only the northeast part of NGC 6888, 



detected p 



( iChu et ah 



asma with a temperature of T pa 2 — 3 MK and possible nitrogen enrichment 
2006) . Sim ilar soft X-ray spect ra were detected for another WBB, S308, with 



ROSAT flWriggdll999f l and X MM- Newton flChu et al.l 120031 ). Unfortunately, none of these 
data provided us with solid grounds to draw a firm conclusion about the physical mechanism 
responsible for the X-ray emission from WBB. Obviously, the X-ray emitting plasma in the 
hot bubble has a moderate temperature (considerably lower than the one expected from a 
shock with velocity equal to that of the stellar wind) and the reason for this could be that 
the thermal conduction operates efficiently in the hot interior. But, it is necessary to have 
solid observational arguments that support or rule out such a theoretical expectation. 

This was the basic motivation for the recent Suzaku observations of NGC 6888 which 
covered the entire object and provided us with the X-ray spectra having good photon statis- 
tics. We report here the results from our analysis of these data. The paper is organized as 
follows. We give some basic information about the WBB NGC 6888 in Section [2J In Section 
|3j we briefly review the Suzaku observations. In Section HI we present results from the global 
spectral models. In Section [51 we discuss the origin of the hot gas in NGC 6888. In Section 
|6l we discuss our results and we present our conclusions in Section 



The Wind-Blown Bubble NGC 6888 



The nebula NGC 6888 around the Wolf-Rayet star WR 136 is the prototype of a broad 
class of ring nebulae around stars with strong stellar w inds. The size o f the optical nebula 
is 12' x 18' and its expansion velocit y is 75 km s -1 ( Chu et al.l 119831 ) . The distance to 
NGC 6888 (WR 136) is d = 1.26 kpc ( Ivan der HuchtlhoOll ). thusTthe dynamical age of the 



object is 30,000 years. The observed X-ray flux in the ROSATPSPC band is 1 2 x 10~ 12 erg 



cm 2 s 1 with an X-ray absorption column density of Nh = 3.2 x 10 21 cm 2 (IWrigge et al. 



19941 ). A large r flux value of 1.86 x 10 12 erg cm 2 s" 



is obtained from the ACS A observations 



of this object (IWrigge et al.l 120051 ). 



The central star, WR 136, is a WN6 star (Ivan der Huchtl 120011 ) with stellar wind 



1600 km s" 1 dPrinia et al.l 11990( 1 and M = 1 5 x 10~ 5 M ? yr 



parameters V W i n d 

For the mass-loss rate, we used the result for WR 136 derived by lAbbott et al.l ( 119861 ). 
logM = —4.5 M yr -1 for d = 1.8 kpc and V W i n d = 2000 km s" 1 , and corrected it for 
the values of the wind velocity and distance adopted here. We made use of the general 
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dependence M oc Vw i ndd 1,5 for the mass-loss values derived from the radio observations 



flPanagia fc Fellil Il975l : IWright fc Barlowlll975r i 



3. Observations and Data Reduction 

NGC 6888 was observed with Suzaku in two runs: June 4 and November 16, 2009, 
centered on the North-East and South- West parts of this object thus covering the entire 
optical nebula. The basic data for our analysis are those from the three active X-ray Imaging 
Spectrometers (XIS) on-board Suzaku. We note that one of them (XIS1) is back-illuminated 
while the other two (XISO, XIS3) are front-illuminated CCDs which allows the extracted 
spectra from units XISO and XIS3 be combined for further analysis 0. We extracted X-ray 
spectra of NGC 6888 from the filtered and screened pipeline events files using the 3x3 and 
5x5 observational modes @ which resulted in 88.8 ksec and 77.3 ksec effective exposures for 
the June and November observations, respectively. The total number of source counts in the 
(0.3 - 4 keV) energy band are 5403 (XIS1) and 4838 (XIS0+XIS3) for the North-East half; 
4330 (XIS1) and 3413 (XIS0+XIS3) for the South- West half of NGC 6888. For the current 
analysis, we extracted one total and two subregion spectra for each of the observed parts of 
NGC 6888. The corresponding X-ray images as well as source and background regions for 
the spectral extraction are shown in Fig. [TJ The circles in red denote the regions around 
bright point sources (No. 5, 6 and 7 in Table 4 of Wrigge et al. 1994) that were excluded 
from the spectral extraction. We note that source No . 5 is in fact the centr al star (WR 136) 



in the nebula. Since it is not a strong X-ray source (jSkinner et al.ll2010l ). the presence of 
a near by point source (No. 6) and the large PSF of Suzaku telescopes prevents extraction 
of valuable information on its X-ray emission. Despite the large 'pixel' of Suzaku images, 
Figure [T] shows that the X-rays of NGC 6888 originate from region inside the optical nebula 
and their spatial distribution is not uniform (e.g., indicating clump like structures). 

We adopted the current calibration database for Suzaku, XIS (20100123) and XRT 
(20080709), to generate the corresponding response matrices and ancillary respo nse files 



We m ade use of standard as well as custom models in version 11.3.2 of XSPEC (lArnaud 



19961 ) for the spectral analysis in this study. 

Finally, to facilitate comparisons of the new results on NGC 6888 with those from the 



1 see The Suzaku Data Reduction Guide: http://heasarc.gsfc.nasa.gov/docs/suzaku/analysis/abc/ 



2 In the 3x3 and 5x5 editing modes, all the pulse heights of the corresponding 9 and 25 pixels centered 
at the event center are sent to the telemetry. The data for each mode come from different times which 
effectively increases the exposure time. See |http: / /heasarc.gsfc.nasa.gov/docs/astroe/prop_tools / suzaku_td| 
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previous X-ray observations of this object, we mention a few basic characteristics of the 
Suzaku observations and similarly for those of ROSAT and ASCA. The Suzaku telescopes 
have spatial resolution of ~ 2 arcmin (expressed as Half-Power Diameter, or HPD). The 
XIS detectors operate in the 0.2 — 12 keV energy range. At energies below 1 keV the 
back-illuminated CCD (XIS1) has much better sensitivity (at least a factor of 2) than the 
front-illuminated CCDs (XISO, 3), while their effective area is about the same at higher 
energies (see § 3 in The Suzaku Technical Description, footnote |2]). On the other hand|f|, the 
energy bandpass and spatial resolution (HPD) were: 0.5 — 12 keV and 2.9 arcmin for ASCA; 
0.12 - 2.4 keV and ~ 30 arcsec for ROSAT. 



4. Global Spectral Models 

4.1. Soft X-ray Spectra 

A general property of the Suzaku spectra of NGC 6888 is that they are rather soft and 
no apparent line features are detected at photon energies > 1.5 keV. Hence, we started our 
analysis with the X-ray spectra in the (0.3 - 1.5 keV) range (re-binned to have a minimum 
100 counts per bin) where most of the X-ray counts (~ 90%) are found. We adopted the 
following approach. 

We used canonical discrete-temperature plasma models to derive the global properties 
of the X-ray emitting plasma in the hot bubble of NGC 6888. Since the ionization history 
of the plasma in the hot bubble is very complex and not well constrained, we considered 
two limiting cases: models that assume plasma in collisional ionization equilibrium (CIE) or 
plasma with non-equilibrium ionization (NIE). For consistency in the atomic data, we made 
use of emission models vapec and vpshock in XSPEC. 

As a first step, we fitted the total spectra of the North-East (combined region 1 and 2 
in Fig. [[]) and South- West (combined region 3 and 4 in Fig. [1]) halves simultaneously (four 
spectra in total). Although the quality of the Suzaku spectra is much bett er compared to that 



of the previous stu dies of the entire object, using ROSAT and ASCA f jWrigge et al.lll994 



Wrigge et aJj|2005| ), we do not find it sufficient to allow us investigate possible variations 



of the X-ray absorption or abundances. Thus, we assumed uniform abundances for the hot 



3 For the ASCA Technical Description see 
http://heasarc.gsfc.nasa.gov/docs/asca/asca_nra06/appendix_e/appendix_e.html 

For the ROSAT Technical Description see |http://heasarc.gsfc.nasa.gov/docs/rosat/ruh/handbook/| and 
|http://heas arc. gsfc.nasa.gov/docs/journal/psp c2. html 
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plasma in NGC 6888 and its X-ray emission is subject to the same interstellar absorption. 
These were free pa rameters in our fits . As a basic set of abundances, we adopted the 
solar abundances of lAnders fc Grevessd (119891 ). We varied only the N, O, Ne, Mg and Fe 
abundances in the spectral fits since these are the elements that provide most of the X-rays in 
the (0.3 - 1.5 keV) energy range. We also explored a different set of (non-solar) abundances 
( Ivan der Hucht et al.lll986l ; see below) that is typical for the stellar wind of the central WN6 
star, since this star is the energy source for the hot bubble in NGC 6888. 

The results from the global model fits to the entire X-ray spectrum of NGC 6888, 
assuming uniform plasma temperature, are given in Table [1] and Fig. [2j We note the good 
quality of the fits and the fact that these fits confirm the early finding; that the X-ray plasma 



in NGC 6888 is relatively 'cool' (IBochkarevI 11988 : IWrigge et all ll994J : Wrigge et all 12005 



Chu et al.l 120061 ). The value of the X-ray absorption is slightly smaller than the one found 



from the analysis o f the entire X-ray spectrum using RO SAT and ASCA (IWrigge et al.lll994 



Wrigge et al.ll2Q05h while the observed flux in the Suazku spectra of NGC 6888 is by 10-20% 



higher. A possible reason for this could be the quality of the spectra and we note that the 
total number of counts of the Suzaku spectra was at least one order of magnitude larger 
than that in the RO SAT and ASCA spectra of NGC 6888. Thus, we believe that the results 
from the current analysis are better constrained. For comparison, Figure [2] also shows the 
background spectra after subtracting the non X-ray background, NXB (the NXB spectra 
were constructed following the procedure in the Suzaku Data Reduction Guide, see footnote 
d in § ED- We see that the source and background spectra are very different (e.g., no N 
VII emission in the latter). This makes us confident about the derived results, despite the 
large 'pixel' of the Suzaku telescopes and the close location of the source and background 
subtraction regions. 

As a next step, we searched for temperature variations of the X-ray plasma in the hot 
bubble of NGC 6888. For this goal, we used IT vpshock model to fit the individual spectra of 
region 1 and 2 in the North-East and of region 3 and 4 in the South- West halves (see Fig. [I]) 
simultaneously (eight spectra in total). All the spectra shared the same abundances and 
X-ray absorption. We note that the derived X-ray absorption (N# = 2.40^ 38 x 10 21 cm -2 ) 
and abundances are consistent (being within la errors) with their values from the spectral 
model fits with uniform plasma temperature and there is no appreciable variation in the 
plasma temperature within the hot bubble: kT= O^Iq;^ (region 1); 0.431q;}q (region 2); 



0.43 + 

o!o9 (region 3); 0.371° 

m (region 4); that is all the values are within the la errors from 



each other. The same is valid for the observed fluxes in the (0.3 - 1.5 keV) energy range: 
Fx = (4.44; 4.69; 4.79 and 4.70) x 10~ 13 ergs cm -2 s _1 for region 1, 2, 3 and 4, respectively. 
The observed spectra overlaid with the best-fit IT shock model {x 2 /dof = 196/232) are 
shown in Fig. [3j Although the X-ray emission of NGC 6888 originates mostly from a non- 
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unifo rm distribution of hot plasma (filaments, clumps; see lWrigge et al.lll994[ IWrigge et al. 
20051 ; see also Fig. [TJ, the properties of the X-ray emitting gas (temperature, flux) are quite 
uniform through the entire hot bubble. So, the high signal-to-noise Suzaku data confirm 
similar findings already reported f r om the analysis of the RO 'SAT and ASCA observations 
of NGC 6888 JWrigge et all ll994J : IWrigge et aDl2005h . But, we caution that the almost 
equal values for the observed flux from regions 1, 2, 3 and 4 do not mean that the plasma 
emissivity in the hot bubble is uniform. As seen from Fig. [TJ regions 1 and 4 comprise much 
smaller part from the nebula than regions 2 and 3 do. Therefore, the plasma emissivity (or 
emission measure) of the X-ray emitting plasma is higher in regions 1 and 4, likely due to 
higher plasma density or higher concentration of clumps and filaments. 

It is interesting to note that although we used simple plasma models to fit the observed 
X-ray spectra of NGC 6888, the derived abundances are fairly consistent between different 
models, and they are generally inconsistent with the typical abundances of a WN star (see 
below). In addition to the two basic models (2T vapec and vpshock) used to fit the entire 
spectrum of NGC 6888 , we ran some more complex ones: (i) a model with a distribution of 
thermal CIE plasma (we adopted our custom model which is similar to c6pvmkl in XSPEC 
but uses the apec collisional plasma for the X-ray spectrum at a given temperature); (ii) a 
model with a distribution of adiabatic shocks with NIE effects taken into account (we made 
use of our custom model for XSPEC which was s uccessfully used in the analysis of the X-ray 
spectra of SNR 1987A - e.g., IZhekov et al.ll2009l and the references therein). We will discuss 
the result from these more complex models below (see § I4.2p . and we only mention that the 
models with a distribution of CIE plasma or NIE shocks show two peaks in the distribution 
of their emission measure with temperature: one for a 'cool' plasma ( kT < 0.4 keV) that 
dominates the X-ray emission from NGC 6888 and another for a 'hot' plasma (kT > 2 keV). 
The values for the plasma temperature of the 2T vapec model are consistent with the values 
of the temperature peaks in the distribution of the CIE plasma. 

For comparison, we a lso ran a vpshock model with a set of typical WN abundances 
(Ivan der Hucht et al.lll986l ). We note that there was no good fit to the observed spectra un- 
less we varied the abundances of N, O, Ne, Mg, and Fe. Figure 0] presents the derived abun- 
dances from various model fits to the entire X-ray spectrum of NGC 6888. All abundances 
are given with respect to the total number density of hydrogen and helium (H+He) which 
facilitates the comparison between cases with completely different basic set of abundances 
(e.g. solar vs. non-solar). On the one hand, this result illustrates the uncertainties of the 
derived abundances if using different plasma models and, on the other hand, it conclusively 
demonstrates that the chemical composition of the X-ray emitting plasma in NGC 6888 is 
quite different from that of the stellar wind of the central star in this object. 



-8 - 



4.2. Very Hot Plasma in NGC 6888? 

In order to search for indications of X-ray emission from hot plasma (kT > 2 keV), we 
expanded the photon range of interest and considered the X-ray spectra of NGC 6888 in the 
(0.3 - 4.0 keV) energy range. This increased the total number of source photons by ~ 10% 
compared to the (0.3 - 1.5 keV) spectra analyzed above. To improve the photon statistics, 
we re-binned the spectra to have a minimum of 200 counts per bin. 

We used two models to fit the spectra: (i) one with a distribution of CIE plasma; (ii) 
and a model with a distribution of NEI shocks (as described in § !4.ip . The quality of the fits 
was very good: x 2 /dof = 122/169 (CIE plasma) and x 2 /dof = 120/168 (NEI shocks). The 
derived values for the X-ray absorption and chemical abundances were entirely consistent 
with the values discussed in our analysis of the soft (0.3 - 1.5 keV) spectra. Figure [5] shows 
the NGC 6888 spectra overlaid with the best-fit model for a distribution of CIE plasma as 
well as the derived distribution of emission measure for the CIE plasma and NEI shocks. 

We see that both models show a bi-modal distribution of the X-ray emitting plasma 
with peaks quite distinguished in the temperature space. The 'cool' plasma (kT < 1 keV) is 
the dominant source of X-ray emission in the (0.3 - 1.5 keV) energy range but provides only 
between (6 - 19)% of the flux at higher energies (1.5 - 4.0 keV), depending on the model used 
(e.g., CIE plasma or NIE shocks). It is thus quite likely that some hot plasma is present 
in the X-ray emitting region of NGC 6888. This plasma is needed to explain the weak but 
statistically significant emission (signal-to-noise ~ 6 — 15) in the (1.5 - 4.0 keV) energy range 
which may have been missed in the ROSAT and ASCA spectra of the entire WBB due to 
their considerably poorer photon statistics, and almost no sensitivity at energies > 2 keV 
for ROSAT. 



5. The Origin of the Hot Gas in NGC 6888 



The two most important results from the model fits to the X-ray spectra of NGC 6888 that 
could be directly related to the origin of the hot gas in this object are: (i) the X-ray plasma 
in NGC 6888 has a relatively low temperature; (ii) the hot bubble abundances are not 
consistent with the ones of the stellar wind of the central star in NGC 6888 (see § H|). 

The latter is justified in two ways as seen from Fig. HI The X-ray abundances are smaller 
than the corresponding canonical WN values (typical for the central star) by (0.5 - 1.5) dex 
on the average. And, the N/O ratio derived from the analysis of the Suzaku data has a 
characteristic value of 2.0 - 3.0 which is much lower than the value of ~ 21.5 for the N/O 



ratio in the WN abundances set (e.g., Ivan der Hucht et al.lll986l ). We note that low values 
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(~ 2) for the N/O ratio were derived from analysis of the optical emission of NGC 6888. 
Also, the N, O a nd Ne abundances derived from X-rays are very similar to those reported 
from the optical (jKwitterlll98ll ; iMoore et al.ll2000t see also Fig. HJ). 



It is thus natural to conclude that the plasma which dominates the X-ray emission from 
NGC 6888 likely originates from the cold optical nebula. Subsequently, this cold plasma has 
been heated to high temperatures through a mechanism that is yet to be identified. This 
heating mechanism must be capable of explaining the relatively low plasma temperature in 
the hot bubble of NGC 6888 (kT < 0.5 keV; see §HJ. 



sec 



As proposed by the stand ard WBB model (1 Weaver et al.lll977l ; for numerical modeling 



Zhekov &: Myasnikovl 11998). the electron thermal conduction is such a physical mecha- 



nism. Alternatively, the models of so called mass- loa ded astronomical flows also suggest a 



relati vely low temperature in the WBB objects (e.g., lHartquist et al. 1 11986c Arthur et al. 
19931 ) but we note that it is not well established which is the actual physical mechanism 
that is responsible for the 'instantaneous' mixing of the hot bubble gas with the mass- loaded 
flow. In either case some important questions must find their answers in order to build a 
more complete global picture of NGC 6888. 

Namely, what is the fate of the shocked stellar wind? Given the age of this WBB (§ [2]), 
we note that the massive WN wind can supply up to one solar mass or so of hot gas with 
distinct non-solar abundances. How can we explain that we find no indications for such a 
plasma in the X-ray spectrum of NGC 6888? And, how can the presence of really high- 
temperature plasma (kT > 2 keV) fit into a general physical picture for WBB? Below, we 
will discuss these issues in some extent. 



6. Discussion 

6.1. Wind-Blown Bubble Hydrodynamics 

In order to have a realistic view on the fate of the shocked WN stellar wind, we need 
to carry a direct comparison between the theoretical X-ray emission and the observed spec- 
tra of NGC 6888. For this goal, we adopted the following procedure. First, we used a ID 
numerical hydrodynamic model to derive the distribution of plasma characteristics (tem- 
perature, density) in the hot bubble. In turn, we constructed the hot plasma distribution 
of emission measure with temperature. Second, we developed a custom model for XSPEC 
that uses thus derived distribution of emission measure to model the X-ray emission from 
NGC 6888. Finally, we used this model to fit the Suzaku spectra and check the consistency 
of the physical model with observations: e.g., whether the model provides correct observed 
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flux, abundances, X-ray absorption. 



We briefly review som e details about the WBB hy drodynamics (see lWeaver et al.lll977 



and numerical modeling in lZhekov &: Myasnikovlll998l ). When a massive stellar wind inter- 



acts with the circumstellar matter, CSM, a two-slab structure will form, bounded by two 
shocks and a contact discontinuity separating the shocked CSM gas from the shocked stellar 
wind. During the evolution of this structure, its outer part will collapse due to the energy 
losses by emission and will appear as optical nebula. On the other hand, its inner part is 
not subject to significant energy losses due to its high temperature and low density but its 
structure will depend on the efficiency of the electron thermal conduction. In the case the 
thermal conduction is suppressed (e.g., due to presence of magnetic field), the hot bubble 
will remain adiabatic and it will contain only the very hot gas of the shocked stellar wind. 
Alternatively, the efficient thermal conduction will heat up part of the cold optical nebula 
which will expand inward (the contact discontinuity will move towards the inner shock). As 
a result, the plasma temperature in the hot bubble will decrease and its mass will increase. 
In this case the hot bubble will consist of two parts both filled with hot gas but with different 
chemical composition. 

To derive the physical parameters of the hot bubble in NGC 6888, we used the ID 
hydrodynamic code of A.V. M yasnikov that was develop ed for modeling the physics of a 



standard WBB (for details see IZhekov fc Myasnikovl Il998l ). This code allows modeling not 



only of adiabatic and radiative WBB but also of the case with efficient thermal conduction. 
It correctly handles interacting gas flows with different chemical compositions. 

We adopted the following basic set of parameters that determine the physics of a WBB: 
(i) Mwn = 1-5 x 10~ 5 M yr -1 and Vwn = 1600 km s" 1 for the stellar wind of the central star 
in NGC 6888 (see §EJ; (ii) M sw = 1.5 x 10~ 4 M yr" 1 and V sw = 10 km s" 1 for the slow wind 
that has been emitted in the previous evolutionary stage of the central star. The interaction 
of such winds results in a WBB with a radius of ~ 2.5 pc at age of 30,000 years. We note 



that the size of the optical nebula NGC 6888 is 12' x 18' flChu et all 11983 : also see Fig. Q}. 
Thus, a spherically-symmetric (ID) bubble with a radius of ~ 2.5 pc has the same volume 
as a prolate ellipsoid with major and minor axes equal to those of the optical nebula for the 
adopted distance to this object (d = 1.26 kpc, §|2j). With this approximation, we calculate 
the distribution of emission measure of the hot bubble and the corresponding normalization 
parameter for the model spectrum in XSPEC (norm = 10" 14 J n e nndV / And 2 ). Then, the 
actual spectral fit to the X-ray spectra can tell us if the theoretical emission measure gives 
the correct value for the observed flux from NGC 6888. 



We fitted the total spectra of NGC 6888 with the theoretical X-ray emission for the 
case of conductive WBB having the basic set of physical parameters. In all model fits using 
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results from our hydro dynamic simulations, the abundances of the sho cked stellar wind had 



values fixed to the typical WN abundances (Ivan der Hucht et al.lll986l ) and only abundances 



of the 'evaporated' gas from the cold shell were allowed to vary. We note that using the basic 
set of parameters provided acceptable fit to the shape of the observed spectra but it grossly 
overestimated the observed flux (by a factor of ~ 50). To reconcile this flux discrepancy, 
it was necessary to reduce the amount of X-ray emitting plasma in the hot bubble. We 
ran a series of hydrodynamic simulations by simultaneously reducing the mass-loss rates of 
the central star in NGC 6888 and of the previously emitted slow wind which results in no 
change of the global geometry of the shock structure. In all these simulations, the wind 
velocities had their values as in the basic set of physical parameters (see above). The case 
with a reduction factor of ~ 4 for the mass-loss rates provided the correct observed flux 
and the derived abundances were a factor of (3 - 5) higher than the values in Fig. |H The 
corresponding spectral fit (x 2 /dof = 266/238) to the Suzaku spectra of NGC 6888 is shown 
in Fig. [6j 

It is worth noting that although the required reduction for the mass-loss rates might 
seem considerable, we find additional evidence that this could well be the case. If we assume 
that the thermal conduction is not efficient and some other mechanism is responsible for 
the relatively low plasma temperature deduced from the analysis of the X-ray emission 
from NGC 6888 (§ @|, the hot bubble in this object will then be purely adiabatic. We 
ran a series of ID hydrodynamic models for this case starting with the same set of basic 
physical parameters as described above. The result was that the shocked stellar WN wind 
plasma had too high a temperature and the models provided strong X-ray flux at energies 
> 1 keV considerably exceeding the one observed. For graphical clarity, we show in Fig. [6] 
the theoretical spectrum of purely adiabatic WBB with mass-loss rates reduced by a factor 
of 2. Note that the normalization (or the flux) of the theoretical spectrum scales with the 
mass- loss rate as cx M 2 . 

Therefore, it is conclusive that independently from which is the heating mechanism for 
the X-ray emitting plasma in NGC 6888, we find indications that the mass-loss rate of the 
central star should be smaller than the value of M^n — 1-5 x 10~ 5 M Q yr _1 adopted here. 
Direct comparisons of the model predictions and the X-ray observations suggests a reduction 
factor of ~ 3 — 4. 



6.2. General Physical Picture 



A general physical picture of the X-ray emission from a WBB must provide a reasonable 
explanation about the origin of the X-ray emitting gas and the mechanism that heats it up 
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to the observed temperatures. Thus, we recall the basic properties of the X-rays from 
NGC 6888. 

The previous observations of the entire WBB with ROSAT and ASCA found that 



the X-ray emission of NGC 6888 comes from clump-like structures (IWrigge et al.l 11994 



Wrigge et al.l 120051 ) . This is also the case as seen in the high signal-to- noise Suzaku data 
(Fig. [p. Our analysis of the Suzaku spectra of NGC 6888 showed that most of the X-rays 
come from relatively cool plasma (kT < 0.5 keV) but there is also some very hot plasma (kT 
> 2 keV; Fig. that is responsible for the observed emission in the (1.5 - 4.0 keV) energy 
region. 

The physical picture that emerges from the analysis of the X-ray data of NGC 6888 is 
then the following. Most of the soft X-rays come from clumps distributed in the hot bubble 
of this object. As indicated by the similar chemical composition of the X-ray emitting plasma 
and the optical shell (e.g., Fig. HJ), these clumps likely originate from the cold but ionized 
optical nebula. Most probably they were formed through vario us dynamic instabilities the 



optical nebula was subject to during its evolution ( e.g. see iGarcia-Segura et al.l Il996a 



Garcia-Segura et al. 1996b ; Stickland fc Stevensl 19981; and multiple sho ck waves may have 



contributed in the case of conductive WBB. IZhekov &: Myasnikovlll998l ). In addition, there 
is some very hot plasma that can be associated with the shocked WN wind of the central 
star. The details of this general physical picture depend on the heating mechanism that 
operates in the hot bubble of NGC 6888. 

Conductive WBB. If the electron thermal conduction is efficient, the energy of the 
shocked stellar wind will be dissipated, part of the outer shell will 'evaporate' and this will 
form a smooth component of relatively hot plasma filling in the inner part of the nebula. 
This case is well represented by the hydrodynamic model of WBB with thermal conduction 
and this model gave an acceptable fit (x 2 /d°f — 266/238) to the observed spectra (see 
§ 16.11 and Fig. [6]). In this picture, the dense ionized clumps that were formed via dynamic 
instabilities are engulfed by the hot bubble and heated up to X-ray temperatures thanks 
to the efficient thermal conduction. In fact, when we added an extra component (using 
model vapec in XSPEC) to the spectral fit based on the hydrodynamic simulations, the 
spectral model matched better the observed spectra (x 2 /d°f — 232/238) and the derived 
abundances had values consistent with those in Fig. HI The derived plasma temperature 
for this additional component (the clumps) was kT = 0.11 — 0.12 keV. The main problem 
in the case of conductive WBB is that there cannot exist in the hot bubble plasma with 
temperatures as high as those derived from the fits to the Suzaku spectra in the (0.3 - 4.0 
keV) energy range (§ !4.2l and Fig. [5]). As a result, practically no photons are produced in the 
(1.5 - 4.0 keV) energy band contrary to what is observed (§ |4.2p . We recall that asymmetric 
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thermal conduction (IZhekov fe Myasnikovl 120001 ) might be a remedy to this problem if a 



global ordered magnetic field is present in NGC 6888. Thus, there might be a small part in 
the hot bubble with a very hot plasma while the rest of it will be filled in with much cooler 
plasma. Unfortunately, the number of detected photons in the (1.5 - 4.0 keV) energy band 
did not allow us to study their spatial distribution and much deeper observations are needed 
for that purpose. 

Adiabatic WBB. If the electron thermal conduction is suppressed due to the presence 
of magnetic field, the hot bubble will be adiabatic and the dense clumps that originate from 
the cold outer shell and penetrate the bubble interior are heated up to X-ray temperatures 
through some mechanism different from thermal conduction. We explored quantitatively 
this case by considering a spectral model consisting of two components: one representing the 
emission from the adiabatic hot bubble (from our hydrodynamic simulations) and another 
one (using model vapec in XSPEC) representing the X-ray emission of the dense clumps. 
We obtained a good fit (x 2 /dof = 183/237) to the Suzaku spectra of NGC 6888 in the 
case of reduced by a factor of 3 mass-loss rates. The temperature of the dense clumps was 
kT = 0.11 — 0.12 keV and the derived abundances were consistent with those in Fig. 0] Then, 
what is the mechanism that heats up the dense clumps? Given the chemical composition 
of the clumps (or the cold nebula), shocks with velocity ~ 300 — 350 km s _1 can provide 
the required temperature (kT = 0.11 — 0.12 keV). But, simple considerations show that the 
rarefied plasma in the hot bubble does not have the necessary power to run such shocks into a 
gas with much higher density than the hot plasma itself (e.g., given the stellar wind velocity 
of 1600 km s _1 , the expected relative velocity between the hot plasma and the clumps is 
maximum of the order of the velocity of the shock). Could it be that the dense clumps 
mix locally with the hot plasma and since the electron thermal conduction is inefficient, the 
energy exchange occurs directly between the heavy particles (ions) of the two components? 
Consequently, the electrons in the mixture are heated as well. Also, the relative velocity 
between the hot plasma and the clumps will favor such a diffusion of the heavy particles. We 
note that the adiabatic WBB does not have the high-temperature caveat of the conductive 
WBB mentioned above. 

But, how can we distinguish between the two cases, conductive or adiabatic WBB? We 
think that there might be some observational differences resulting from these two heating 
mechanisms. The electron thermal conduction is a very efficient dissipative mechanism, thus 
the heating of the cold gas of the clumps will occur almost 'instantaneously'. On the other 
hand, the ion-ion energy exchange (thanks to diffusion of the hot ions into the cold gas) is 
much slower process. We recall that this different efficiency is due to the much faster energy 
exchange between the light particles (the characteristic time is proportional to the square root 
of the particle's mass; e.g., ISpitzerlll962l ). Thus, we may expect that there will exist clumps 
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in the hot bubble that have temperature intermediate between that of the optical outer shell 
and the X-ray clumps. These clumps will be source of UV emission. Therefore, if spectral 
observations of NGC 6888 detect strong lines in the 1,000-2,000 A spectral range (e.g., of 'hot' 
ions like NV, OVI), this will indicate that the thermal conduction is considerably suppressed 
in this object. Alternatively, if clumps with no strong UV emission are detected, this will 
favor the thermal conduction as the heating mechanism in NGC 6888. Another obsery a tiona l 



possibility also exists, although it is not quite feasible for the moment. iNordon et all (120091 ) 
reported a narrow radiative recombination continuum (RRC) in the X-ray spectrum of the 
planetary nebula BD+60°3639. They interpret this RRC emission feature as result when 
bare C VII ions penetrating from the hot bubble into the cold optical nebula recombine 
with cool electrons. In general, the same can be expected when the hot ions diffuse into the 
cold clumps as proposed above and RRC of various ions can be expected. But, a crucial 
experiment for distinguishing between detection from non-detection of these RRC features 
must wait for future observations with the next generation of X-ray telescopes since such an 
experiment requires high spectral resolution X-ray observations which are not feasible for 
NGC 6888 at the moment. 

All this, as well as the suggested reduction for the mass-loss rate of the central WN star, 
emphasizes the importance to carry out some global analysis of the entire system NGC 6888: 
central star, optical and X-ray nebula (hot bubble). Namely, a self-consistent modeling of the 
optical-UV emission of the central star (WR 136) and the optical nebula and confronting 
it with observations can provide us with a valuable piece of information about the wind 
velocity and mass loss, temperature distribution in the stellar wind, its abundances as well 
as the density and abundances of the optical nebula. This will put much tighter constraints 
on the physical parameters determining the physical conditions in the hot bubble (and its 
X-ray emission). In turn, it will result in a better constrained general physical picture of the 
whole system. A suitable example of adopting such a global approach in th e analysis of the 



optical-UV emission of a complex system is given in iGeorgiev et all (120081 ) for the case of 
the planetary nebula NGC 6543. 



6.3. Ionization History of the Hot Plasma 

The ionization history of the plasma in the hot bubble of a WBB is quite complex. 
Consider the standard WBB model. In this case, the thermal conduction is efficient and we 
have two parts in the hot bubble filled correspondingly with the shocked stellar wind and 
the evaporated gas from the cold optical nebula. During the evolution of the WBB, the 
temperature in the hot bubble decreases with time as well as does its density. Thus, each 
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new parcel of gas, crossing the inner shock or evaporating from the cold nebula, will have 
its own ionization history. This will depend on the time that parcel of gas entered the hot 
bubble, that is on the physical parameters of the hot plasma and their subsequent evolution 
until the moment of observations. We note that it is likely (although it depends on the 
physical parameters for each studied object) that the plasma on both sides in vicinity of the 
contact discontinuity will be almost in CIE while that close to the inner shock or to the outer 
boundary in cold nebula will be in state of NEI. In purely adiabatic WBB, the ionization 
history of the plasma is similar but then the hot bubble contains only the shocked stellar 
wind. This picture becomes even more complicated in objects like NGC 6888 where the 
X-ray emission is not spatially uniform which indicates a presence of clump like structures. 

Unfortunately, it is not feasible to follow the ionization history of the plasma in the hot 
bubble in the theoretical models or to deduce valuable information about it from the X-ray 
observations. The latter is particularly valid for the case of CCD X-ray spectra. For that 
reason, our approach in this study was to consider two limiting cases: plasma models with 
CIE and models that, although in a simplified but well defined manner, take into account the 
NEI effects. Such an approach illustrates better the uncertainties in the derived properties 
(e.g., temperature, abundances, emission measure) of the X-ray emitting plasma. 



7. Conclusions 

In this work, we presented an analysis of the Suzaku data on NGC 6888: the wind-blown 
bubble around the Wolf-Rayet star WR 136. For the first time, X-ray spectra of the entire 
nebula are obtained that have a very good quality allowing a detailed comparison between 
theory and observations. The basic results and conclusions are as follows. 

1. The X-ray spectra of the WBB NGC 6888 are soft and most of the emission is in the 
(0.3 - 1.5 keV) energy range although some emission (~ 10% of the observed flux) 
is found at higher energies (1.5 - 4.0 keV). The spectral fits require a relatively cool 
plasma with kT < 0.5 keV but much hotter plasma with temperature kT > 2.0 keV is 
needed to match the observed hard X-ray emission. We find no appreciable temperature 
variations within the hot bubble of NGC 6888. 

2. One of the important results from the spectral fits is that the derived abundances 
(N, O, Ne) are very much the same as those of the optical nebula. Also , they are 



considerably different from the WN abundances ( Ivan der Hucht et al.lll986l ) assumed 
'typical' for the wind of the central star in the nebula. This indicates a common origin 
of the X-ray emitting gas and the outer cold shell, that is most of the X-ray plasma 
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(likely concentrated in clumps as the X-ray images indicate) has flown into the hot 
bubble from the optical nebula. 

3. A direct comparison (in XSPEC) between ID hydrodynamic models and the X-ray 
spectra of NGC 6888 suggests a reduced mass-loss rate (~ 4 — 5 x 10~ 6 M Q yr _1 ) of 
the central star in order to provide the correct value of the observed flux. We note that 
this figure is appreciably low for Wolf-Rayet stars and thus emphasizes the need of a 
global modeling of the entire system: central star, optical nebula and the hot bubble. 

4. The general physical picture that emerges from the current analysis is the following. 
Most of the X-rays come from clumps distributed in the hot bubble of NGC 6888. 
These clumps originate from the cold optical nebula and they formed through various 
dynamic instabilities during the evolution of the wind-blown bubble. The shocked WN 
wind of the central star is likely the source of the very hot plasma (kT > 2 keV) in this 
wind-blown bubble. If the electron thermal conduction is efficient, this can naturally 
explain the relatively low plasma temperature of most of the X-ray emitting plasma. 
Alternatively, the hot bubble in NGC 6888 will be adiabatic and the cold clumps are 
heated up to X-ray temperatures likely by energy exchange between the heavy particles 
(e.g., a local diffusion of the rarefied hot plasma into the dense clumps). 
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Fig. 1.— NGC 6888 images: optical (Palo mar Observatory Sky Survey; upper panel); 
Suzaku XIS1 (0.3 — 4.0 keV) for the June 2009 (lower left panel) and November 2009 (lower 
right panel) observations. R.A. (J2000) and decl. (J2000) are on the horizontal and vertical 
axes, respectively. The boxes in green denote the extraction regions for the source (numbered 
from 1 through 4) and background spectra (labeled 'Bkg'). The circles in red mark the regions 
around point sources that were excluded from the spectral extraction. The circled plus sign 
gives the position of the central star WR 136 in NGC 6888 . In lower panels, the circle 
in green (labeled 'HPD') illustrates the half-power diameter (~ 2 arcmin) of the Suzaku 
telescopes. The X-ray images are re-binned (4x4 original pixels) and wavelet smoothed for 
presentation. For clarity, overlaid are some contours (in black) of the optical image. Note 
that the lower panels show the raw X-ray images, thus, their brightness scale is different due 
to the different non X-ray background. 
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Fig. 2. — NGC 6888 background-subtracted spectra overlaid with the best-fit IT shock model 
(Tabled]). In each upper panel, the XIS1 spectrum is in black and the XISO+3 spectrum is 
in red (lower curve). The North-East and South- West halves of the bubble were observed 
correspondingly in June and November 2009. For comparison, the lower panels show the 
corresponding background spectra after subtracting the non X-ray background. 
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Fig. 3. — NGC 6888 background-subtracted spectra from the four extraction regions in Fig. [T] 
overlaid with the best-fit IT shock model. In each panel, the XIS1 spectrum is in black and 
the XISO+3 spectrum is in red (lower curve). 
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Fig. 4. — The abundance values of N, 0, Ne, Mg and Fe derived from the global model 
fits to the entire X-ray spectrum of NGC 6888 . The abundance values are with respect to 
the total number density of H+He and in units of 12 + logX (X is any of N, 0, Ne, Mg 
and Fe). The symbols with error bars denote results from the fits with vpshock and solar 
abundances - squares; vpshock and WN abundances - asterisks; distribution of CIE plasma 
- diamonds (the values from the 2T vapec model fits are not given since they are within 
the plotted error bars); distribution of NEI shocks - crosses. Errors are la valu e s from the 
fits. The abundanc es values derived from analysis of the optical data by iKwitterl ( )198ll ) and 



Moore et al.l (120001 ) are shown w ith filled and empty circle s, respectively. For reference, the 



typical abundances of WN stars ( Ivan der Hucht et al.lll986l ) are given by empty squares with 
horizontal bars (these bars have no meaning and are used only for graphical clarity). 
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Fig. 5. — The (0.3 - 4.0 keV) background-subtracted spectra of NGC 6888 (re-binned to 
have a minimum of 200 counts per bin) overlaid with the best-fit model with a distribution 
of CIE plasma. In each panel, the XIS1 spectrum is in black and the XISO+3 spectrum is in 
red (lower curve). The derived distribution of emission measure (normalized to have a total 
sum of unity) is shown in the third panel: that of the CIE plasma in black (solid line) and 
the one of the NEI shocks in blue (dashed line). 
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Fig. 6. — NGC 6888 background-subtracted spectra as in Fig. [2] overlaid with the best-fit 
model from the ID hydrodynamic simulations for conductive WBB with mass losses reduced 
by a factor of ~ 4. In each panel, the XIS1 spectrum is in black and the XISO+3 spectrum 
is in red (lower curve). For comparison, the model spectra from a purely adiabatic WBB are 
shown in green (XIS1 spectrum) and blue (XISO+3 spectrum). See text for details. 
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Table 1. Global Spectral Model Results 



Parameter 2T vapec vpshock 



X 2 /dof 

(10 21 cm- 2 ) 
kTi (keV) 
kT 2 (keV) 
EMi (10 55 cm" 3 ) 
EM 2 (10 55 cm" 3 ) 
Abundances 
r (10 10 cm" 3 s) 
F x (10~ 12 ergs cm -2 s _1 ) 



172/232 


177/235 


1 Q2+ - 42 


2 07 +0 - 47 

z - u ' -0.46 


o 1 5 +0 - 01 

U ' iO -0.01 


0.48±°; 22 






80.8™ 




1 o +0 - 17 

L - u -0.14 




Varied 


Varied 



2.05 (18.4) 1.92 (17.4) 



Note. — Results from simultaneous fits to the XIS1 
and XISO+3 spectra of NGC 6888 from the Suzaku ob- 
servations in June and Nov 2009. Tabulated quanti- 
ties are the neutral hydrogen absorption column den- 
sity (N#), plasma temperature (kT), emission mea- 
sure (EM = J n e riHdV), shock ionization time-scale 
(r = n e t), the X-ray flux (F x ) in the 0.3 - 1.5 keV 
range followed in parentheses by the unabsorbed value. 
Only abundances of N, O, Ne, Mg and Fe were varied 
in the fits (see Fig. HJ). Errors are the la values from 
the fits. 



